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Abstract 

Silicon detectors in particle physics experiments at the new accelerators or 
in space missions for physics goals will be exposed to extreme radiation con- 
ditions. The principal obstacles to long-term operation in these environments 
are the changes in detector parameters, consequence of the modifications in 
material properties after irradiation. 

The phenomenological model developed in the present paper is able to ex- 
plain quantitatively, without free parameters, the production of primary de- 
fects in silicon after particle irradiation and their evolution toward equilibrium, 
for a large range of generation rates of primary defects. Vacancy-interstitial 
annihilation, interstitial migration to sinks, divacancy and vacancy-impurity 
complex (VP, VO, V2O, dOi and CiC s ) formation are taken into account. 
The effects of different initial impurity concentrations of phosphorus, oxygen 
and carbon, as well as of irradiation conditions are systematically studied. 
The correlation between the rate of defect production, the temperature and 
the time evolution of defect concentrations is also investigated. 

PACS: 

29.40 Pe: Semiconductor detectors 

61.80 Az: Theory and models of radiation defects 

61.70 At: Defect formation and annealing 

Keywords: silicon, detectors, radiation damage, kinetics of defects, an- 
nealing processes 

1 Introduction 

The use of silicon detectors in high radiation environments, as to be expected 
in future high energy accelerators or in space missions, poses severe problems due 
to changes in the properties of the material, and consequently influences the per- 
formances of detectors. 

As a consequence of the degradation to radiation of the semiconductor material, 
an increase of the reverse current due the reduction of the minority carrier lifetime, 
a reduction of the charge collection efficiency and a modification of the effective 
doping, due to the generation of trapping centres, are observed in the detector 
characteristics. 
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In this paper, the effects of irradiation conditions and various initial impurities 
in the starting material are discussed in the frame of the phcnomcnological model 
able to explain quantitatively defect production and evolution toward stable defects 
during and after irradiation in silicon. The model supposes three steps. 

In the first step, the incident particle, having kinetic energies with values in the 
intermediate up to high energy range, interacts with the semiconductor material. 
The peculiarities of the interaction mechanisms are explicitly considered for each 
kinetic energy. 

In the second step, the recoil nuclei resulting from these interactions lose their 
energy in the lattice. Their energy partition between displacements and ionisation 
is considered in accord with the Lindhard theory (Q and authors' contributions 
§)■ 

A point defect in a crystal is an entity that causes an interruption in the lat- 
tice periodicity. In this paper, the terminology and definitions in agreement with 
M. Lannoo and J. Bourgoin || are used in relation to defects. We denote the 
displacement defects, vacancies and interstitials, as primary point defects, prior to 
any further rearrangement. After this step the concentration of primary defects is 
calculated. 

The mechanisms of interaction of the incident particle with the semiconductor 
lattice, accompanied by displacement defect production have been discussed in some 
papers, see, e.g. references ||, ||, |(|. The incident particle produces, as a consequence 
of its interaction with ions of the semiconductor lattice, cascades of displacements. 

In silicon, vacancies and interstitials are essentially unstable and interact via 
migration, recombination, and annihilation or produce other defects. 

The concentration of primary defects represents the starting point for the follow- 
ing step of the model, the consideration of the annealing processes, treated in the 
frame of the chemical rate theory A review of previous works about the problem of 
the annealing of radiation induced defects in silicon can be found, e.g. in Reference 

I- 

Without free parameters, the model is able to predict the absolute values of the 
concentrations of defects and their time evolution toward stable defects, starting 
from the primary incident particle characterised by type and kinetic energy 

The first two steps have been treated extensively in previous papers (J, ||, [|, 
where the concentration of primary defects has been calculated. In this paper, the 
third step is discussed extensively and represents a generalisation of the previous 
results published in references [Q, || including also the carbon contributions to defect 
kinetics. 

The influence and the effects of different initial impurity concentrations of phos- 
phorus, oxygen and carbon as well as of the irradiation conditions were system- 
atically studied. Vacancy-interstitial annihilation, interstitial migration to sinks, 
vacancy - impurity complexes (VP, V2, VO, V2O, CiOi, CiC s ) - only the stable de- 
fects confirmed experimentally in silicon for high energy physics applications were 
considered. The correlation between the rate of defect productions, the temper- 
ature and the time evolution of the defect concentrations was also investigated. 
Some conclusions about the possibilities to obtain semiconductor materials harder 
to radiation are given. 

2 Production of primary defects 

The basic assumption of the present model is that the primary defects, vacancies 
and interstitials, are produced in equal quantities and are uniformly distributed in 
the material bulk. They are produced by the incoming particle, as a consequence 
of the subsequent collisions of the primary recoil in the lattice, or thermally (only 
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Prenkel pairs arc considered). The generation term G is the sum of two components: 

G = G R + G T (1) 

where Gr accounts for the generation by irradiation, and Gt for thermal generation. 
The concentration of the primary radiation induced defects per unit fluence CPD 
in silicon has been calculated as the sum of the concentration of defects resulting 
from all interaction processes, and all characteristic mechanisms corresponding to 
each interaction process, using the explicit formula (see details, e.g. in references 
1,0): 

CPD (E) = ^i- fy(^r) L{E Rl ) Sl dn = J-^l^NIEL (E) (2) 
1 ' 2E Si J Z-< \d£l) it3i V Rl ' Sl N A 2E Si y ' K> 

where E is the kinetic energy of the incident particle, N$i is the atomic density 
in silicon, A$i is the silicon atomic number, E$i - the average threshold energy 
for displacements in the semiconductor, Em - the recoil energy of the residual 
nucleus produced in interaction i, L(ERi) - the Lindhard factor that describes the 
partition of the recoil energy between ionisation and displacements and (der/dO), 

- the differential cross section of the interaction between the incident particle and 
the nucleus of the lattice for the process or mechanism i, responsible in defect 
production. Na is Avogadro's number. The formula gives also the relation with the 
non ionising energy loss (NIEL) the rate of energy loss by displacement (dE/dx) n i 

10- 

In Figure 1, the kinetic energy dependence of CPD for different particles is 
presented: for pions the values are calculated in accord with equation (1) and 
arc from reference [fL2[ , and for the other particles these are evaluated from the 
published NIEL calculations, as follows: for protons - reference || for neutrons 

- reference jlJJ , for electrons - reference || and for photons from || . The main 
source of errors in the calculated concentration of primary defects comes from the 
modelling of the particle - nucleus interaction and from the number and quality of 
the experimental data available for these processes. 

Due to the important weight of annealing processes, as well as to their very 
short time scale, CPD is not a measurable physical quantity. In silicon, vacancies 
and interstitials are essentially unstable and interact via migration, recombination, 
and annihilation or produce other defects. In the simplifying hypothesis of random 
distribution of CPD for all particles, used in the present paper, the identity of 
the particle is lost after the primary interaction, and two different particles could 
produce the same generation rate (Gr) for vacancy-interstitial pairs if the following 
condition: 

G R = [(CPD) part a {Ex)] ■ % art .a{El) = [(CPD)^ (E 2 )\ ■ % art . 2{ E 2) (3) 

is fulfilled. Here, $ is the flux of particles (a) and (b) respectively, and E\ and E 2 
their corresponding kinetic energies. 

3 The kinetics of radiation induced defects 

Silicon used in high energy physics detectors is n-typc high resistivity (4-5-6 
YSl cm) phosphorus doped FZ material. 

The effect of oxygen in irradiated silicon has been a subject of intensive studies 
in remote past. In the last decade a lot of studies have been performed to investigate 
the influence of different impurities, especially oxygen and carbon, as possible ways 
to enhance the radiation hardness of silicon for detectors in the future generation 
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Figure 1: Energy dependence of the concentration of primary defects on unit fluence 
induced by protons, pions, electrons, photons and neutrons in silicon - see text for details. 



of experiments in high energy physics - see, e.g. references 15, [T(|. Some people 
consider that these impurities added to the silicon bulk modify the formation of 
electrically active defects, thus controlling the macroscopic device parameters. Em- 
pirically, it is considered that if the silicon is enriched in oxygen, the capture of 
radiation-generated vacancies is favoured by the production of the pseudo-acceptor 
complex vacancy-oxygen. Interstitial oxygen acts as a sink for vacancies, thus re- 
ducing the probability of formation of the divacancy related complexes, associated 
with deeper levels inside the gap. 

The concentrations of interstitial oxygen Oi and substitutional carbon C, in sil- 
icon are strongly dependent on the growth technique. In high purity Float Zone 
Si, oxygen interstitial concentrations are around 10 15 cm -3 , while in Czochralski Si 
these concentrations can reach values as high as 10 18 cm -3 . Because Czochralski 
silicon is not available in detector grade quality, an oxygenation technique devel- 
oped at BNL produces Diffusion Oxygenated Float Zone in silicon, obtaining a 
Oi concentration of the order 5xl0 17 cm -3 . These materials can be enriched in 
substitutional carbon up to 1.8xl0 16 cm~ 3 . 

After the irradiation of silicon, the following stable defects have been identified 



(see References § |17|): Si t , VP, VO, V 2 , V 2 0, C.O,, C u C l C s . The pre-existing 
thermal defects and those produced by irradiation, as well as the impurities, are 
assumed to be randomly distributed in the solid. An important part of the va- 
cancies and intcrstitials annihilate. The sample contains certain concentrations of 
impurities, which can trap intcrstitials and vacancies respectively, and form stable 
defects. 

Vacancy-interstitial annihilation, interstitial migration to sinks, divacancy, va- 
cancy and interstitial impurity complex formation are considered. The role of phos- 
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phorus, oxygen and carbon is taken into account, and the following stable defects : 
VP, VO, V2, V2O, CiOi, Ci, CiC s are considered. Other possible defects as V3O, 



V2O2, V3O3 1 18 , are not included in the present model. 

The following picture describes in terms of chemical reactions the mechanisms 
of production and evolution of the defects considered in the present paper: 

V + I ^annihilation (4) 
/ — > sinks (5) 



VO is the A centre. 



VP is the E centre. 



V + O^VO (6) 

V + P p VP (7) 



V + V ^ V 2 (8) 

V + VO ^? V 2 (9) 

/ + C s ^ d (10) 

d + Oi^dOi (11) 

A + I^O (12) 

C t + C s ^ CiC s (13) 

Some considerations about the determination of the reaction constants are given in 
references ||] . 

The reaction constants Ki (i = 1, 3 -j- 9) have the general form: 

tfi =Cuexp(-Ei/k B T) (14) 

with ^ the vibration frequency of the lattice, Ei the associated activation energy 
and C a numerical constant that accounts for the symmetry of the defect in the 
lattice. 

The reaction constant related to the migration of interstitials to sinks could be 
expressed as: 

K 2 = av\ 2 exp {-E 2 /k B T) (15) 

with a the sink concentration and A the jump distance. 

The system of coupled differential equations corresponding to the reaction scheme 
(4) -7- (13) cannot be solved analytically and a numerical procedure was used. 

The following values of the parameters have been used: E\=E 2 = 0.4 eV, E 3 
= 0.8 cV, E 4 = 1.4 cV, E 5 = 1.1 eV, E 6 = 1.3 eV, E 7 = 1.6 cV, E 8 = 0.8 cV, E 9 
= 1.7 eV, v = 10 13 Hz, A = 10 15 cm 2 , a = 10 10 cm~ 2 , in accord with standard 
constants characterising the silicon material, see for example the books of Lannoo 
and Bourgoin || or Damask and Dienes @. 
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4 Results, discussion, comparison with experimen- 
tal data and predictions 



The formation and time evolution of stable defects depends on various factors, 
e.g. the concentrations of impurities pre-existent in the sample, the rate of genera- 
tion, and the temperature. 

In Figures 2 a -4- f , and 3 a 4- f the formation and time evolution of the vacancy- 
oxygen, vacancy-phosphorus, divacancy, divacancy-oxygen, carbon interstitial - oxy- 
gen interstitial and carbon interstitial.- carbon substitutional are modelled in silicon 
containing different initial concentrations of phosphorus, oxygen and carbon, and 
for two very different rates of generation. 




Figure 2: Time dependence of the concentrations of a)VO, b)VP, c)V2, d)V20, e)dOi 
and f) CiC s induced in silicon with the following concentrations of impurities: (1): 10 14 
P/cm 3 , 2xl0 15 O/cm 3 , and 3xl0 15 C/cm 3 ; (2): 10 14 P/cm 3 , 10 16 O/cm 3 , and 3xl0 15 
C/cm 3 ; (3):10 14 P/cm 3 , 5xl0 16 O/cm 3 , and 310 15 C/cm 3 ; (4): 10 14 P/cm 3 , 4xl0 17 O/cm 3 , 
and 3xl0 15 C/cm 3 ; (5)): 10 15 P/cm 3 , 10 16 O/cm 3 , and 3xl0 15 C/cm 3 (6): 10 14 P/cm 3 , 
10 16 O/cm 3 , and 3xl0 4 6 C/cm 3 by a generation rate G R = 7x10 s VI pairs/cm 3 s, at 20°C. 

In Figure 2, the evolution of defect concentrations during high rate irradiation 
(Gr = 7xl0 8 VI pairs/cm 3 s) is presented. This corresponds, in the model hypoth- 
esis, to order of magnitude of the radiation levels estimated for the forward tracker 
region at the future LHC accelerator. 
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The increase of the initial oxygen concentration in silicon, from 2xl0 15 O/cm 3 to 
4x1 17 O/cm 3 , conduces, after ten years of operation in the field characterised by a 
high and constant generation rate, to the increase of the concentrations of VO and 
CiOi centres, and to the decrease of the concentrations of V2, VP and CiC s ones. 
With the increase of oxygen concentration, a variation of the V2O generation rate is 
observed, so that, from the studied cases, the maximum concentration for this defect 
is obtained for 5xl0 16 O/cm 3 initial oxygen. The increase of initial phosphorus is 
seen in the increase of concentration of VP centres, while the increase of initial 
carbon concentration has important consequences on the concentrations of CiC s 
centres. It is interesting to observe that in almost all cases, an equilibrium in 
reached between generation and annealing, and a plateau is obtained in the time 
dependence of the concentrations. The slowest is, in this respect, V2O, that has the 
highest binding energy. 

As underlined above, vacancy-oxygen formation in oxygen enriched silicon is 
favoured in respect to the generation of V2, V2O and VP. This is an important 
feature that could be used for detector applications, determining the decrease of 
the leakage current J8|. At high oxygen concentrations, the concentration of VO 
centres attains a plateau during the 10 years period considered. 
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Figure 3: Same as Figure 2, Gr = 200 VI pairs/cm 3 s. 
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The other extreme situation corresponds to a generation rate Gr = 200 VI 
pairs/cm 3 s, equivalent with a rate of defect production by protons from the cosmic 
ray spectra in the orbit near the Earth, at about 400 Km, as will be the position 
of the International Space Station. The same concentrations of P, O and C have 
been considered as pre-cxistent in silicon as in Figure 2. For this generation rate, 
the increase of the oxygen concentration produces the decrease of the concentration 
of all centres, with the exception of the VO concentration, that, at these rates, it 
is not influenced by the oxygen content, and of the CiOi concentration, where an 
increase is observed. As could be seen from Figure 3, as a consequence of the small 
rate of generation rate of vacancy - interstitial pairs, after ten years of operation the 
equilibrium between generation and annealing is not reached, the concentrations of 
defects being, with the exception of VP (that has a relatively low binding energy) , 
slightly increasing functions of time. 

All curves, both from Fig. 2 and Fig. 3, have been calculated for 20°C temper- 
ature. Thermal generation has been taken into account in both cases, although it 
is important only for the silicon exposed to low rates of defect production. 

The influence of the generation rate of primary defects on the concentration of 
stable defects and on their time evolution has also been investigated. In Figure 4 
a -T- f , the time evolution of the VO, VP, V2, V2O, CiOi and CiC s concentrations 
is presented for six decades of generation rates of defects ((7r=7x10 5 4- 7xl0 10 VI 
pairs/cm 3 s), for silicon containing the following concentrations of impurities: 10 14 
P/cm 3 , 10 16 O/cm 3 , and 310 15 C/cm 3 at 20°C temperature. At small times, the 
curves corresponding to different generation rates are all parallel and equidistant 
in a log-log representation. Starting from the highest generation rates, they start 
to increase slower (V2, VP), attain a plateau V2O, or event start to decrease (VO, 
CiOi, CiC s . The maximum attained can be the same, independent on the gener- 
ation rate as is the case of VO concentration, or could be generation dependent 
CiOi, CiC s . 

The temperature is another important factor determining the time evolution of 
defects. In Figure 5 a f f and 6 a f f, the effect of the temperature is studied 
for the same generation rates of primary defects as in Figures 2 and 3 respectively, 
for five temperatures: 20°C, 10°C, 0°C, -10°C, and -20°C. The concentrations of 
pre-existing impurities in silicon are as follows: 10 14 P/cm 3 , 10 16 O/cm 3 , and 
3xl0 15 C/cm 3 . The decrease of the temperature decreases the generation rate 
of stable defects, with the exception of VP, where the highest values correspond 
to the lowest temperature. The maximum values of the VO concentration are 
temperature independent; the values of the concentration on the plateau decrease 
with the decrease of temperature for CiOi and CiC s and increase for VzO. 

For silicon exposed to low rates of defect production the same amount of time 
(Figure 6), the process of defect production in slowed down. The most unexpected 
time dependence is for V2, that has the highest values at the lowest temperature. 

In a previous paper [gj, we demonstrated in concrete cases the importance of 
the sequence of irradiation process, considering that the same total fluence can be 
attained in different situations: the ideal case of instantaneous irradiation, irradia- 
tion in a single pulse followed by relaxation, and respectively continuos irradiation 
process. As expected, after instantaneous irradiation the concentrations of defects 
are higher in respect with "gradual" irradiation. 

After this analysis, the specific importance of the irradiation and annealing 
history (initial material parameters, type of irradiation particles, energetic source 
spectra, flux, irradiation temperature, measurement temperature, temperature and 
time between irradiation and measurement) on defect evolution must be to under- 
line. 

The model predictions have been compared with experimental measurements. 
A difficulty in this comparison is the insufficient information in published papers 
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Figure 4: Time dependence of the concentrations of a) VO, b)VP, c)V2, d)V20, e)dOi 
and f) dCs induced in silicon with: 10 14 P/cm 3 , 10 16 O/cm 3 , and 3xl0 15 C/cm 3 by 
continuous irradiation. 



regarding the characterisation of silicon, and on the irradiation parameters and 
conditions for most of the data. 

It was underlined in the literature |2(J that the ratio of VO to VP centres in 
electron irradiated silicon is proportional to the ratio between the concentrations of 
oxygen and phosphorus in the sample. For electron irradiation, in Ref. [|TJ a linear 
dependence of the Vi versus VO centre concentration has been found experimen- 
tally. In the present paper, the ratio of concentrations of V2 to VO centres and VO 
to VP ones has been calculated in the frame of the model, for the material with 
the characteristics specified in Ref. j2^], and irradiated with 12 MeV electrons, in 
the conditions of the above mentioned article. The time dependence of these two 
ratios is represented in Figure 7. Annealing is considered both during and after 
irradiation. It could be seen that for the ratio of V2 and VO concentrations the 
curves corresponding to different irradiation fluences are parallel, while the ratio of 
VO to VP concentrations is fluence independent, in the interval 2xl0 13 ~ 5.5xl0 14 
cm~ 2 , in good agreement with the experimental evidence. The ratio between V2 
and VO concentrations is determined by the generation of primary defects by irra- 
diation, while the ratio between VO and VP concentrations is determined by the 
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Figure 5: Time dependence of the concentrations of a) VO, b)VP, c)V2, d)V20, e)dOi 
and f) dCs induced in silicon with: 10 14 P/cm 3 , 10 16 O/cm 3 , and 3xl0 15 C/cm 3 by 
continuous irradiation, in the same conditions as in Figure 2, at 20°C, 10°C, 0°C, -10°C, 
and -20° C. 



concentrations of oxygen and phosphorus in silicon. 

Our estimations are also in agreement with the measurements presented in ref- 
erence ^2|, after electron irradiation, where defect concentrations are presented as 
a function of the time after irradiation. In Figure 8, both measured and calculated 
dependencies of the VO and VP concentrations are given. The irradiation was 
performed with 2.5 MeV electrons, up to a fiuence of 3xl0 16 c/cm 2 . A good agree- 
ment can be observed for these concentrations. The dependencies put in evidence 
the important role played by the carbon-related defects. The relative values are 
imposed by the arbitrary units of experimental data. 

Also, a good agreement has been obtained between the absolute values of con- 
centrations of VP + V2 and CiC s predicted by the model, and the experimental 
results after neutron irradiation at a total fiuence of 5.67xl0 13 cm -2 , reported in 
reference ^3|. The calculated 1.5xl0 13 cm -3 and 4.1xl0 12 cm~ 3 concentrations for 
VP + V2 and CiC s respectively, are in accord with the values of l.lxlO 13 cm~ 3 
and 3.8xl0 12 cm -3 , measured experimentally. For the VO concentration, a poorer 
concordance has been obtained. 
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Figure 6: Time dependence of the concentrations of a) VO, b)VP, c)V2, d)V20, e)dOi 
and f) dCs induced in silicon with: 10 14 P/cm 3 , 10 16 O/cm 3 , and 3xl0 15 C/cm 3 by 
continuous irradiation, in the same conditions as in Figure 3, at 20°C, 10°C, 0°C, -10°C, 
and -20° C. 



5 Summary 

A phcnomenological model that describes silicon degradation due to irradiation 
from the point of view of the kinetics of produced defects toward equilibrium was 
developed. 

The production of primary defects (vacancies and interstitials) in the silicon 
bulk was considered in the frame of the Lindhard theory, and the peculiarities of 
the particle - silicon nuclei interaction were taken into account. 

The mechanisms of formation of stable defects and their evolution toward equi- 
librium was modelled, and the concentrations of defects were calculated solving nu- 
merically the system of coupled differential equations for these processes. Vacancy- 
interstitial annihilation, interstitial migration to sinks, vacancy-impurities com- 
plexes (VP, VO, V%0, CiOi, CiC s ), and divacancy formation were considered in 
different irradiation conditions, for different concentrations of impurities in the ini- 
tial semiconductor, in the temperature range —20 -j- +20°C. 

The calculated results suggest the importance of the conditions of irradiation, 
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a) 




Figure 7: Time dependence for a) V2/VO and b) VO/VP concentrations calculated 
for silicon with 1.4xl0 14 P/cm 3 , 5xl0 17 O/cm 3 , and 3xl0 15 C/cm 3 by 12 MeV electron 
irradiation, with the flux 5.8xl0 10 e/cm 2 s, up to the fluences: 2xl0 13 e/cm 2 , 4xl0 13 e/cm 2 , 
8xl0 13 e/cm 2 , 2xl0 14 e/cm 2 and 5.5xl0 14 e/cm 2 , followed by relaxation (see reference ]2l| . 



temperature and annealing history. 

The model supports the experimental studies performed to investigate the influ- 
ence of oxygen in the enhancement of the radiation hardness of silicon for detectors. 
The VO defects in oxygen enriched silicon are favoured in respect to the other stable 
defects, so, for detector applications it is expected that the leakage current decreases 
after irradiation. At high oxygen concentrations, this defect saturates starting from 
low fluences at high generation rates of defects. 
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